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, Abstract. The stellar initial mass function and the stellar lifetimes are basic ingredients of chemical evolution 

T— I ■ models, for which different recipes can be found in the literature. In this paper, we quantify the effects on chemical 

^ ' evolution studies of the uncertainties in these two parameters. We concentrate on chemical evolution models for 

. the Milky Way, because of the large number of good observational constraints. Such chemical evolution models 

have already ruled out significant temporal variations for the stellar initial mass function in our own Galaxy, 
with the exception perhaps of the very early phases of its evolution. Therefore, here we assume a Galactic initial 
mass function constant in time. Through an accurate comparison of model predictions for the Milky Way with 
' carefully selected data sets, it is shown that specific prescriptions for the initial mass function in particular mass 

, ranges should be rejected. As far as the stellar lifetimes are concerned, the major differences among existing 

' prescriptions are found in the range of very low-mass stars. Because of this, the model predictions widely differ 

i";^ ' for those elements which are produced mostly by very long-lived objects, as for instance ^He and ^Li. However, it 

l^i is concluded that model predictions of several important observed quantities, constraining the plausible Galactic 

Q ' formation scenarios, are fairly robust with respect to changes in both the stellar mass spectrum and the stellar 

^ I lifetimes. For instance, the metallicity distribution of low-mass stars is nearly unaffected by these changes, since 

C/3 . its shape is dictated mostly by the time scale for thin-disk formation. 

■ 
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1. Introduction Pettini 2001; Centurion et al. 2003; Dietrich et al. 2003a, 

b; Prochaska, Howk & Wolfe 2003; Tolstoy et al. 2003; 

The formation and evolution of galaxies is one of the out- Dessauges-Zavadsky et al. 2004; D'Odorico et al. 2004). 

standing problems of astrophysics. In the last decade, a this framework, galactic chemical evolution models can 

great deal of observational work has shed light on the regarded as useful tools to discriminate among dif- 

production and distribution of chemical elements inside fgrent scenarios of galaxy formation. In fact, as stressed 

the Galaxy (e.g. Edvardsson et al. 1993; Cayrel 1996; j^any times in the literature, abundances and abundance 

Nissen & Schuster 1997; Gratton et al. 2000; Chen et al. ratios play a major role as cosmic clocks and give hints 

2003; Gratton et al. 2003; Ivans et al. 2003; Reddy et the time scales of structure formation and evolution 

al. 2003; Zoccah et al. 2003; Akerman et al. 2004), often (wheeler, Sneden & Truran 1989; Matteucci & Frangois 

leading to an evolutionary scenario much more compli- 1992- Matteucci 2001) 

cated than assumed in many models. Even more recently, order to build up a chemical evolution model, it 

abundance data have accumulated for external galaxies at jg necessary to define the initial conditions and the ba- 

both low and high redshift, thus providing precious in- physical laws governing the evolution of the system 

formation on the chemical evolution of different types of ^^^^^^ ^he whole galactic lifetime. In short, one needs to 

galaxies and on the early stages of galaxy evolution (e.g. gpg^ify whether the system is closed or open (should any 

inflow/outflow of gas occur), the chemical composition of 

Send offprint requests to: D. Romano the gas from which the computation starts and that of 
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any infalling material, the stellar birthrate function and 
the stellar evolution and nucleosynthesis. In particular, 
the stellar birthrate is expressed as the product of two in- 
dependent functions, the star formation rate (SFR) and 
the stellar initial mass function (IMF). The first is gener- 
ally expressed as a function of time only, while the second, 
which describes the stellar mass distribution at birth, is 
likely to be universal (Kroupa 2002; but see e.g. Jeffries 
et al. 2004) and not vary as a function of time (Chiappini, 
Matteucci & Padoan 2000). 

The free parameters that one introduces in the model 
simply reflect our poor understanding of the basic phys- 
ical processes governing the formation and evolution of 
galactic structures. However, having a number of obser- 
vational constraints formally larger than that of the free 
parameters allows us to restrict the range of variation of 
the parameters themselves and gain useful insight into the 
mechanisms of galaxy formation and evolution (Matteucci 
2001). This is the case for our own Galaxy and will become 
soon the standard for an increasing number of external 
galaxies, thanks to the capabilities of modern telescopes 
and instrumentation. 



In an epoch where the uncertainties in the data have 
become really small, it is worth trying to consider 'the- 
oretical error bars' as well. An attempt to do so has re- 
cently been done by Romano et al. (2003), who compare 
the evolution of light elements predicted by two indepen- 
dent models of chemical evolution for the Milky Way and 
try to ascertain the origin of the differences in the model 
predictions. In the present paper we intend to assess the 
uncertainties in the model predictions which arise when 
exploring different prescriptions for the IMF and the stel- 
lar lifetimes in a model for the Milky Way. To this pur- 
pose, we adopt a chemical evolution model which has been 
proven to successfully reproduce the main observational 
features of the solar neighbourhood (Chiappini, Matteucci 
& Gratton 1997) and change the assumptions on the stel- 
lar IMF. As a result, we set up a range of possible varia- 
tions for several predicted quantities. Then, we repeat the 
same kind of analysis by changing the prescriptions for 
the stellar lifetimes. 



The paper is organized as follows. In Section 2, we re- 
view the general features of the IMF which have emerged 
over the years and emphasize differences and similari- 
ties among the various IMFs adopted in this work. In 
Section 3, we discuss the prescriptions for the stellar life- 
times. In Section 4, we describe the chemical evolution 
model for the solar vicinity. In Section 5 we present model 
results. Finally, Section 6 is devoted to a critical discussion 
of the problem and some conclusions are drawn. Notice 
that the paper is structured in such a way that it is easy 
to concentrate on only one among the proposed topics, 
while skipping the others, if one wishes to do so. 



2. The stellar initial mass function 

2.1. The adopted parametrizations 

The most widely used functional form for the IMF is an 
extension of that proposed by Salpeter (1955) to the whole 
stellar mass range: 

<;^Salpeter(m) = ^Salpeter m"(^+'^\ 

where x = 1.35 and 



J m 



?7l</'Salpeter(m)dm = 1. 



Here the IMF is by number; miow = 0.1 Mq, mu 



100 



Mq and Salpeter — 0.17. The above formula equivalently 

reads 



<<?Salpeter(m)dm = 1 



if using the IMF by mass, (/?Saipeter(w) cx m~^'^^. In what 
follows we always use the IMF by mass. 

More realistic, multi-slope expressions better describe 
the luminosity function of main sequence stars in the solar 
neighbourhood, that is what is actually observed (Tinsley 
1980; Scalo 1986; Kroupa, Tout & Gilmore 1993; Scalo 
1998; see the original publications for details): 

-1.0 



I'Tinsley ' 



VTinsley("^) = { -^Tinsley W 
^Tinsley m' 



-1.3 
-2.3 



if TO < 2 M0 
if 2 < tu/Mq < 10 
if TO > 10 Mq, 



■<*Tmsley — 0.21, <^Tinsley — 0.26, ^Tinsley — 2.6; 



</'Scalo86('7l) 



•2^Scalo86 m' 

^Scaio86 rn' 



-1.35 
-1.70 



if TO < 2 Mq 

if TO > 2Mq, 



Scalo 86 



0.19, 



5Scalo86 



0.24 



(notice that we adopt a simplified two-slope approxima- 
tion to the actual Scalo 1986 formula, similarly to what is 
done in Matteucci & Frangois 1989); 



VKroupa(m) 



Kroupa ^ 
^Kroupa ^ 



-0.3 
-1.2 



, TO 



-1.7 



if TO < 0.5 Mq 

if 0.5 < to/Mq < 1 
if TO > IMq, 



-^Kroupa 



¥'Scalo98(w) 



0.58, ^Kroupa = 'i 

•a?scaio98 rn~ 

■^Scalo98 m' 



Kroupa — 0.31; 
0-2 if TO < 1 Mp 



Scalo 98 m' 



-1.7 
-1.3 



© 



if 1< to/Mq < 10 

if TO > 10 Mq, 



^Scalo98 = ^Scalo98 - 0.39, '^Scalo98 - 0.16. AU of them 

are considered in the present work. The normalization is 
always performed in the mass range 0.1-100 Mq. 

More recently, a lognormal form has been suggested 
for the low- mass part of the IMF (to < 1 Mq), eventually 
extending into the substellar regime (Chabrier 2003): 



¥'Chabrier(?n) = < 



=^^^habrier 6 



'Chabrier W 



if TO < 1 MfT 



-(1.3±0.3) 







if TO > 1 Mq 



D. Romano et al.: Quantifying the uncertainties of diemicai evolution studies 3 

Table 1. I(mi, TO2) (see text for a definition of this quantity) for different mass intervals for the various IMFs 
considered in this study. 
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Fig. 1. a) Stellar IMF according to Salpeter (1955; triangles), Tinsley (1980; stars), Scalo (1986; full circles), Kroupa 
et al. (1993; squares), Scalo (1998; empty circles) and Chabrier (2003; dots). Here (p{m) is the IMF by mass and 
(p{m) oc {x = 1.35 in the case of Salpeter's IMF), except for Chabrier (2003), where a lognormal form for the 

low-mass domain [m < 1 M0) is suggested instead (see text). Here the Chabrier IMF also has x = 1.3 as the exponent 
in the m > 1 Mq mass domain, b) Same as panel a), but with ip{m) divided by the corresponding Salpeter value 
for each given mass. This allows a first sight comparison of the various mass distributions expected according to the 
different IMFs with respect to the 'standard' Salpeter choice. 



According to Chabrier (2003), the IMF depends weakly on 
the environment, except perhaps for early star formation 
conditions. Values of rric = 0.079 Mq and a = 0.69 well 
characterize the IMF for single objects belonging to the 
Milky Way disk. For m > 1 Mq, we assume a power-law 
exponent x equal to either 1.3 or 1.7. In the latter case, 
we study both an IMF truncated at mjow = 0.1 Mq and 
one extending down to toiow — 0.001 Mq, i.e., into the 
brown dwarf (BD) domain. The normalization constants 

are: ^/chabrier — 0.85, ^Chabrier — 0.24 for X = 1.3, TOlow 
= 0.1 Mq; M^habrier ^ 1-16, .^Chabrier ^ 0.32 for X = 1.7, 
miow = 0.1 Mq; M^habrier ^ 1-06, ^Chabrier ^ 0.30 for X 

= 1.7, miow = 0.001 Mq. 



The question of whether the IMF is universal or varies 
with varying star-forming conditions has been long de- 
bated in the past. Combining recent IMF estimates for 
different populations in which individual stars have been 
resolved unveils an extraordinary uniformity of the IMF 
(Kroupa 2002), although some room for exceptions is left 
(e.g., Aloisi, Tosi & Greggio 1999). Explaining the chemo- 
photometric properties of elliptical galaxies has often re- 
quired an IMF biased towards massive stars (e.g. Arimoto 
& Yoshii 1987). However, the actual IMF slope cannot be 
inferred from direct observations in this case. 
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Table 2. I(mi, 1712) for different mass ranges in case of Chabrier (2003) IMF with x = 1.7. The normalization is 
performed either in the mass range 0.1-100 Mq (second column) or in the mass range 0.001-100 Mq (third column). 
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Fig. 2. Stellar IMF according to Chabrier (2003). A power 
law is assumed for m > 1 Mq, with an exponent x = 1.3 
(dots) or X = 1.7 (pentagons). For m < 1 Mq, a lognor- 
mal form is proposed, eventually extending into the BD 
domain (open pentagons). 



2.2. Differences and similarities among different 
parametrizations 

In Fig. ^ we compare Salpeter (1955; triangles), Tinsley 
(1980; stars), Scalo (1986; full circles), Kroupaet al. (1993; 
squares), Scalo (1998; empty circles) and Chabrier (2003; 
dots) IMFs. For the Chabrier IMF the dots for m > 1 Mg 
display the power-law form with an exponent x = 1.3. For 
all these IMFs, the normalization is performed in the mass 
range 0.1-100 Mq. 

It is immediately seen that a Salpeter or a Scalo (1986) 
IMF predicts many more stars at the very low end of 
the distribution than a Scalo (1998) or a Chabrier IMF. 
Tinsley's IMF lies somewhere in the middle, while it pre- 



dicts the highest fraction of stars in the mass range 2-10 
Mq. These features appear more clearly in Fig.^^j where 
the IMFs are divided by the corresponding Salpeter value 
for each given mass. 

In order to quantify the relative weights of stars be- 
longing to different mass ranges according to different 
IMFs, we compute the following quantities: 

I(toi,TO2)= / (p{m)dm, 

J nil 

for each of the IMFs discussed above. The results are listed 
in Table 1. The integration limits, mi and m2, are properly 
chosen in order to allow a meaningful comparison when 
discussing the role played by stars belonging to different 
mass intervals in polluting the interstellar medium (ISM) 
with different chemical elements while varying the IMF 
prescriptions (Section 5.1). 

In Fig. 13 we show the effect of steepening Chabrier 
(2003) IMF for m > 1 Mq. A value of a: = 1.7 (pentagons) 
is now assumed instead of cc = 1.3 (dots). A steeper IMF 
for stars more massive than a few solar masses seems to be 
more likely, on the grounds of recent results by Kroupa & 
Weidner (2003). Those authors argue that field IMFs for 
early-type stars must be steeper than the Salpeter approx- 
imation, owing to the fact that a Salpeter power law well 
describes the distribution of stellar masses for local star 
formation events in star clusters. A steep field-star IMF 
thus arises naturally because of the power-law cluster mass 
function according to which star clusters distribute. 

Fig. 13 also shows the behaviour of Chabrier IMF for 
m < 1 Mq. It is seen that the mass distribution flattens 
below m ~ 1 Mq, reaches a peak around m ~ 0.1 Mq 
and then decreases smoothly for m < 0.1 Mq. As a conse- 
quence, less than 10% of the stellar mass should be found 
in form of BDs (see also Table 2). 

In Section 5.1 we illustrate the results obtained for the 
solar vicinity when the different parametrizations listed 
above are adopted. 

3. The stellar lifetimes 

Different authors traditionally adopt different prescrip- 
tions for the stellar lifetimes. For instance, Matteucci 
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Fig. 3. a) Stellar lifetimes as a function of the initial mass of the star: Maeder & Meynet (1989; dots); Tinsley (1980; 
pentagons); Schaller et al. (1992; stars); Padovani & Matteucci (1993; empty circles); Kodama (1997; filled circles); 
Portinari et al. (1998; asterisks), b) A zoom of the 1-10 Mq stellar mass range. 



and coworkers have usually assumed stellar lifetimes from 
Maeder & Meynet (1989): 



-[^Q-0.6545 logm+1 
I^Q — 3.71og m+1.35 
-|^Q-2.511ogm+0.77 
2^0-1.78 logm+0.17 
-[^Q-0.861ogm-0.94 

1.2 + 0.003 



for m < 1.3 Mq 
for 1.3 < to/Mq < 3 
for 3 < ui/Mq < 7 
for 7 < ui/Mq < 15 
for 15 < ui/Mq < 60 
for m > 60 Mq, 



with Tm in units of Gyr. Notice that they extrapolated the 
values for m < 1.3 Mq and m > 60 Mq, since Maeder & 
Meynet (1989) do not give any formula for these mass 
ranges. The adopted extrapolation to larger masses is 
compatible with the values reported by Maeder & Meynet 
for stars of 85 and 120 Mq (their table 2). 

These stellar lifetimes differ widely from other formula- 
tions available in the literature, especially in the low-mass 
range (to < 1 Mq). Tinsley (1980) and Tosi (1982 and 
subsequent papers) have: 



> 8.6 
8.6 > T„ > 0.64 
0.64 > > 0.087 
0.087 > T„, 0.0155 
~ 0.0155 > Tm 0.003 



for m < 1 Mq 
for 1 < to/Mq < 2 
for 2 < to/Mq < 4 
for 4 < to/Mq < 10 
for 10 < m/M© < 50 



(see Tinsley 1980 for references). 
Kodama (1997) reports: 

50 for TO < 0.56 Mq 

j^g (0.334-^1.790-0. 2232 X (7. 764-logm))/o. 1116 



1.2m 



-1.85 



for TO < 6.6 Mq 
0.003 for TO > 6.6 A//q, 



similarly to what suggested by Padovani & Matteucci 
(1993 and references therein): 

160 for m < 0.60 Mq 

(0.334-^1.790-0.2232 X (7.764-log m)) /0.III6 

for TO < 6.6 Mq 

1.2 TO-1-^5 + 0.003 for TO > 6.6 Mq. 

It is seen that ancient very low-mass stars (0.6 < to/Mq < 
0.9) die in a Hubble time according to Matteucci & 
Frangois (1989) extrapolation of Maeder & Meynet (1989) 
formulae, whereas they have a longer life according to other 
authors. On the contrary, following Maeder & Meynet 
(1989) stars with 1 < to/A/q < 2 live longer than in any 
other study. 

In Fig. 13 we compare all the above-mentioned 
parametrizations, as well as values obtained from Geneva 
(Schaller et al. 1992) and Padua (Portinari et al. 1998 
and references therein) stellar tracks. In the low- and very 
low-mass stellar mass domain substantial differences are 
found, while smaller, though not negligible, differences ex- 
ist at higher masses. Detailed analyses of the results ob- 
tained when integrating these differences over the Galactic 
lifetime and IMF are reported in Section 5.2. 

4. The chemical evolution model for the solar 
neighbourhood 

In order to follow the chemical evolution of the solar vicin- 
ity, we adopt the two-infall model of Chiappini et al. 
(1997). According to this model, the halo and part of the 
thick-disk population form out of a first infall episode on 
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a short time scale, while the thin disk accumulates much 
more slowly during a second independent infall episode. 
The rate of accretion of matter is 



dSinf(0 

dt 



where Sinf (i) is the mass surface density of the infalling 
primordial matter at time t. The parameters <max — 1 Gyr, 
th = 0.8 Gyr and td = 7 Gyr are the time of maximum 
infall onto the disk and the time scales for mass accre- 
tion onto the halo/thick-disk and thin-disk components, 
respectively. They are fixed by the request of reproducing 
a number of observational constraints (Matteucci 2001). 
The coefficients A and B are derived from the condition 
of reproducing the current total mass surface density at 
the solar position (Chiosi 1980). Obviously, the coefficient 
B must be zero for t < tmax- 
The SFR is 



S(^Gal) 

m 



k-l 



proportional to both the total mass and gas surface 
densities. Here G{t) is the normalized gas surface den- 
sity, Sgas(i)/S(tGai), and E(iGai) is the total mass sur- 
face density at the present time, ^oai — 13.7 Gyr^. A 
gas exponent k — 1.5 guarantees a good agreement be- 
tween model predictions and observations (Chiappini et 
al. 1997). Moreover, it is found to agree with inferences 
from observations (Kennicutt 1998) and A^-body simula- 
tions (Gerritsen & Icke 1997). 

The star formation efficiency, I'it), is set to 2 Gyr^^ 
during the halo/thick-disk phase and to 1.2 Gyr~^ dur- 
ing the thin-disk phase to ensure the best fit to all the 
observed features of the solar vicinity, unless the gas sur- 
face density drops below a critical threshold, Sg^g ~ 7 
Mq pc~^. In this case i/ = and the star formation 
ceases. This naturally explains the existence of a gap in 
the SFR between the halo/thick-disk and the thin-disk 
phase (Gratton et al. 1996, 2000; Fuhrmann 1998, 2004). 
Moreover, it delays the beginning of the star formation in 
the halo to the time at which the critical gas density can 
be reached. 

The instantaneous recycling approximation is relaxed, 
i.e. the stellar lifetimes are taken into account in details. 
Stellar nucleosynthesis is taken from (i) van den Hoek 
& Groenewegen (1997) for low- and intermediate-mass 
stars (their case with metallicity-dependent mass loss effi- 
ciency along the asymptotic giant branch) ; (ii) Woosley & 
Weaver (1995) for massive stars; (iii) Thielemann, Nomoto 
& Hashimoto (1993) for Type la supernovae (SNela); (iv) 
Jose & Hernanz (1998) for classical novae. The carbon 
yields from massive stars in the mass range m > 40 Mq 
are multiplied by a factor of three (arguments are given in 
Chiappini, Romano & Matteucci 2003a). Stellar produc- 
tion and/or destruction of the light elements deuterium, 

^ Notice that here we adopt a younger Galaxy age to allow 
for consistency with the recent WMAP data on the universe 
age (see also Romano et al. 2003). 



■^He and ^He are included in the model according to 
Dearborn, Steigman & Tosi (1996), GaUi et al. (1997) and 
Chiappini, Renda & Matteucci (2002) prescriptions. ^Li 
production/destruction is accounted for following Romano 
et al. (2001, 2003) and references therein. 

The prescriptions on the IMF and the stellar lifetimes 
are changed according to the above discussions (Sections 2 
and 3). 

5. Quantifying the uncertainties in chemical 
evolution model predictions 

5.1. Changing the IMF 

In this section we discuss the results obtained by varying 
the prescriptions on the stellar IMF in the chemical evolu- 
tion code for the solar vicinity. The purpose is to associate 
errors due to uncertainties in the stellar IMF to chemical 
evolution model results. 

During the Galactic lifetime, many successive stellar 
generations form according to the adopted IMF (see e.g. 
Tables 1 and 2). As a result, for each IMF choice the com- 
posite stellar population, determining the global chemical 
properties of the Galaxy, is more or less enriched in stars 
falling in a given mass range. Hence the model predictions 
on specific quantities related to the given mass range vary 
when changing the IMF prescriptions. 

As a first example. Fig. 0] displays the theoretical 
G-dwarf metallicity distribution predicted by the model 
under different assumptions on the IMF and with fixed 
stellar lifetimes (namely those by Maeder & Meynet 
1989). The theoretical distributions are convolved with 
a Gaussian to account for both the observational and 
the intrinsic scatter. The adopted total dispersion is a 
= 0.15 dex in [Fe/H] (Arimoto, Yoshii & Takahara 1992; 
Kotoncva et al. 2002). The various panels show, clockwise 
from top left, the theoretical distributions expected when 
assuming a Scalo (1986), Salpeter (1955), Chabrier (2003) 
or Scalo (1998) IMF (thick lines). In particular, Chabrier's 
IMF has X — 1.7 as the exponent of the power-law for 
TO > 1 Mq. For all these IMFs, the normalization is per- 
formed in the mass range 0.1-100 Mq. The distribution 
obtained assuming Kroupa et al. IMF looks indistinguish- 
able from that obtained adopting Chabrier IMF. Similarly, 
the distribution obtained with Tinsley IMF looks like that 
predicted with Scalo's (1998). Therefore, we do not show 
them. The position of the peak and, in general, the shape 
of the distribution arc rather insensitive to the assumed 
IMF. They rather depend on the adopted time scale of disk 
formation (e.g. Chiappini et al. 1997). In particular, only 
a few stars are found at [Fe/H] < —0.7 dex, independently 
of the choice of the IMF. Here we should point out that 
our model predictions are for thin disk stars, whereas the 
observed metallicity distribution cannot avoid some con- 
tribution from thick-disk stars, especially at low metallic- 
ities. 

Also shown are data from Rocha-Pinto & Maciel (1996; 
thin solid histograms) and J0rgensen (2000; thin dotted 
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Table 3. Type la SN rates and Type la to Type II SN 
rate ratios at the present time predicted by the model 
with different IMFs. All the IMFs are normalized to the 
mass range 0.1-100 Mq. The value of the parameter A, 
determining the stellar mass fraction belonging to binary 
systems with the right characteristics to give rise to SNIa 
explosions, is the same independently of the adopted IMF. 
Observed values are from van den Bergh & Tammann 
(1991) and Cappellaro et al. (1997). 



IMF A 


(century^^) 


Ri./Rii 


S55 0.05 


0.57 


0.34 


T80 0.05 


1.16 


0.30 


S86 0.05 


0.45 


0.46 


K93 0.05 


0.67 


0.38 


S98 0.05 


1.03 


0.26 


COS 0.05 


0.69 


0.34 


Observed 


0.6 h'^ 


0.15-0.27 



S55-Salpeter (1955); T80-Tinsley (1980); S86 - Scalo (1986); 
K93 - Kroupa et al. (1993); S98 - Scalo (1998); C03 - Chabrier 
(2003) with X = 1.7. 

Table 4. Same as Table 3, but changing the value of A 
to recover the same SNIa rate at the present time with 
different IMFs. 



IMF A 


_Ria (century ^) 


Ri./Rn 


S55 0.04 


0.47 


0.27 


T80 0.02 


0.46 


0.17 


S86 0.05 


0.45 


0.46 


K93 0.035 


0.45 


0.25 


S98 0.023 


0.47 


0.12 


C03 0.032 


0.43 


0.23 


Observed 


0.6 h-' 


0.15-0.27 



S55-Salpeter (1955); T80-Tinsley (1980); S86 - Scalo (1986); 
K93 - Kroupa et al. (1993); S98 - Scalo (1998); C03 - Chabrier 
(2003) with X = 1.7. 



histograms). The two distributions look quite different. 
In the case of J0rgensen's, the peak is shifted to higher 
metallicities and the low-metallicity tail is almost absent, 
similarly to what found by Wyse & Gilmore (1995). This 
is because of the different metallicity calibrations and the 
different corrections to the raw data applied by the au- 
thors. Indeed, in Rocha-Pinto & Maciel (1996) the metal- 
licity scale is biased towards metal-poor stars (Martell & 
Laughlin 2002). This determines a shift of the peak of the 
distribution towards lower metallicities. The actual peak 
position is likely to be more near —0.15 dex than —0.25 
dex (H. Rocha-Pinto, private communication), thus bring- 
ing model predictions and observational distributions in 
better agreement. 

It is worth reminding that iron originates mostly from 
SNIa explosions occurring in binary systems with an 
intermediate-mass primary. Only one third of its solar 
abundance is related to SNII events, occurring on much 
shorter time scales. Therefore, the stars responsible for the 



observed behaviour of the iron abundance belong mostly 
to the 1.5-8 Mq stellar mass range, that is the one char- 
acterizing Type la SN progenitors in our model. 

For each stellar generation, the mass fraction in form 
of binaries having the right characteristics to end up as 
SNela is a free parameter of the model, whose value is kept 
constant in time. Following Matteucci & Greggio (1986), 
the rate of Type la SNe is: 



i?ia(t) = A 



0.5 



Jd^ 



druB, 



where (p{mB) is the IMF for the total mass of the binary 
system; /(/i) is the distribution function for the mass frac- 
tion of the secondary component (/i = 7712/7713); is the 
lifetime of the secondary; ms^ — 3 Mq and tobm = 16 
Mq (see Matteucci & Recchi 2001 for a review and al- 
ternative formulations). The parameter A is fixed by the 
request of reproducing the SNIa rate currently observed in 
the disk. Moreover, the ratio i?ia/-Rii between the Type la 
and II SN rates in the solar vicinity should be reproduced. 

In Table 3 we list the present-day values of i?ia and 
Ria/Ru which are predicted by assuming different IMFs, 
while keeping the value of A constant (the often quoted 
value of 0.05 for A is adopted), li h — 0.7 and a factor 
of two of uncertainty in the observations is assumed, it 
turns out that both a Tinsley (1980) and a Scalo (1998) 
IMF provide only a marginal agreement with the observed 
SNIa rate. The reason for this is that, among all the stud- 
ied IMFs, Tinsley's (1980) and Scab's (1998) are those 
allocating the highest mass fraction to the stellar mass 
range 1.5-8 Mq, typical of SNIa progenitors. 

Table 4 illustrates what happens if the value of A is 
adjusted so as to bring the value of i?ia in agreement 
with the observations when changing the IMF. The cor- 
responding G-dwarf metallicity distributions are shown in 
Fig.^(dashed areas). Lowering A from 0.05 to 0.04, 0.023 
and 0.032, in case of Salpeter (1955), Scalo (1998) and 
Chabrier (2003) IMF, respectively, results in narrower, 
more peaked distributions. Nevertheless, after convolution 
with a Gaussian accounting for the observational and in- 
trinsic error, the differences almost completely level off. 
The agreement with the observations is always satisfac- 
tory. The Chabrier IMF also gives a good fit to the Ri^/Rw 
ratio currently observed in the solar vicinity. For the Scalo 
(1986) IMF, a better value for the i?ia/i?ii ratio can be 
obtained by further reducing the A value. The Kroupa et 
al. (1993) IMF is an equally valid choice. However, no firm 
conclusions can be drawn on the grounds of the SN rates, 
given the high uncertainty level still affecting the data. In 
what follows, for each IMF we discuss the results obtained 
with the corresponding A value listed in Table 4, unless 
otherwise specified. 

In Table 5 we display the abundances of ''He, C, N, 
O, Ne, Mg, Si, S, Fe and the global metallicity, Z, pre- 
dicted by the model at the time of the Solar System for- 
mation 4.5 Gyr ago, under different prescriptions on the 
IMF. Theoretical values are compared to observed ones 
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Table 5. Abundances in the Protosolar Cloud as predicted by the model at t = 9.2 Gyr with different assumptions on 
the IMF. The abundances are by number in loge(X), except for helium (Y) and global metallicity {Z), for which the 
abundances by mass are reported. Model predictions are compared to observed photospheric solar abundances, unless 
otherwise stated. 



IMF 


Y 


C 


N 





Ne 


Mg 


Si 


g 


Fc 


Z 


S55a 


0.264 


8.53 


8.33 


9.12 


8.27 


7.55 


7.81 


7.39 


7.68 


0.028 


T80a 


0.273 


8.56 


8.46 


8.98 


8.07 


7.39 


7.82 


7.40 


7.71 


0.024 


S86a 


0.260 


8.36 


8.04 


8.77 


7.84 


7.14 


7.52 


7.11 


7.51 


0.014 


K93a 


0.264 


8.46 


8.18 


8.87 


7.97 


7.27 


7.62 


7.21 


7.55 


0.017 


S98a 


0.271 


8.64 


8.49 


9.18 


8.36 


7.64 


7.87 


7.44 


7.68 


0.033 


C03a 


0.265 


8.48 


8.21 


8.89 


8.00 


7.30 


7.64 


7.23 


7.55 


0.018 


C03b 


0.263 


8.44 


8.15 


8.85 


7.95 


7.25 


7.60 


7.19 


7.54 


0.017 


Observations 


GS98 


0.275 ±0.01 


8.52 ± 0.06 


7.92 ± 0.06 


8.83 ± 0.06 


8.08 ± 0.06- 


7.58 ± 0.05 


7.55 ± 0.05 


7.33 ±0.11 


7.50 ± 0.05 


0.017 


HOI 




8.59 ± 0.11 


7.93 ± 0.11 


8.74 ± 0.08 




7.54 ± 0.06 


7.54 ± 0.05 




7.45 ± 0.08 





APOl 8.69 ± 0.05 

AP02 8.39 ± 0.04 

A04 8.66 ±0.05 7.84 ± 0.06 0.0126 



^ Coronal data. 



S55 - Salpeter (1955); T80 - Tinsley (1980); S86 - Scalo (1986); K93 - Kroupa et al. (1993); S98 - Scalo (1988); C03 - Chabrier (2003) with x = 
1.7. The letter refers to the adopted normalization range: a stands for 0.1-100 Mq; b for 0.001-100 Mq. GS98 - Grevesse & Sauval (1998); HOI - 
Holweger (2001); APOl - AUende Prieto et al. (2001); AP02 - AUende Prieto et al. (2002); A04 - Asplund et al. (2004). 



from Grevesse & Sauval (1998), Holweger (2001), AUende 
Prieto, Lambert & Asplund (2001, 2002) and Asplund 
et al. (2004). One may wonder whether the use of the 
Sun as representative of the chemical composition of the 
ISM of the local disc 4.5 Gyr ago is suitable for com- 
parison with GCE model predictions. Indeed, the empir- 
ical age-metallicity relationship for solar neighbourhood 
stars shows a large dispersion (Edvardsson ct al. 1993), 
with known planet-harbouring stars being systematically 
more metal rich than stars without planets (Ibukiyama & 
Ariinoto 2002, their figure 21). On the other hand, a re- 
cent reappraisal of the chemical composition of the Orion 
nebula suggests heavy element abundances only slightly 
higher (~ 0.1 dex) than the solar ones, in agreement with 
GCE model predictions, challenging the view that the Sun 
has abnormally high metal abundances (Esteban et al. 
2004). 

It is seen that a model adopting Tinsley (1980) or Scalo 
(1998) IMF is in better agreement with the solar helium 
mass fraction (see also Romano et al. 2003), but over- 
estimates the overall metallicity of the gas and the iron 
content. This is due to the high stellar mass fraction dis- 
tributed over the mass range 2-40 Mq according to these 
IMFs (see Table 1), coupled with the high helium and 
global metal yields predicted in this mass range by stellar 
evolution theory (see figures 2.9, 2.10 of Matteucci 2001). 
Trying to reduce the predicted solar iron abundance by 
further lowering the A parameter is not a viable solution 
(although at first glance it could appear as a very simple 
and promising one), since a too low Ria/Ru ratio would be 
obtained in this case, in disagreement with the available 
data. Notice that a better agreement with the "'He data 
can be achieved by model S86a if adopting the ^He yields 
recently computed by Meynet & Maeder (2002), which in- 
clude mass loss and rotation in self-consistent stellar evo- 
lutionary models. Indeed, a value of 0.272 for Y at the 
time of Sun formation is predicted in this case, in very 



good agreement with the observed solar value (Chiappini, 
Matteucci & Meynet 2003b). 

Salpeter's IMF also predicts far too high values for 
the solar abundances of iron, oxygen and metals in gen- 
eral, because of its high percentage of massive stars (see 
also Tosi 1988, 1996). In particular, with this IMF and 
the star formation law adopted here, there happens to 
be no way to reconcile the theoretical solar abundances 
with the observed ones, unless one requires a SF process 
so unefficient that severe discrepancies with other obser- 
vational constraints do arise. In particular, by using v ~ 
0.2 Gyr~^ in the SFR expression reported in Section 4, 
the model matches the actual metallicity of the Sun, but 
underestimates the present-day stellar mass density in the 
solar vicinity. In fact, in this case the expected stellar mass 
density at the present time turns out to be 5]stars(iGai) — 
22 Mq pc~^, to be compared with an observed value of 
35 ± 5 Mq pc-2. 

The nitrogen solar abundance is overestimated by all 
the models, except by the model assuming Scalo's (1986) 
IMF. This same model also predicts a solar carbon abun- 
dance in good agreement with the observations, if one be- 
lieves that the solar photospheric abundance from Allendc 
Prieto et al. (2002) gives the best estimate of carbon in 
the Sun. A model assuming Kroupa et al. or Chabrier's 
IMF also agrees with the CNO data, but only at the 2-a 
level. However, it should be noticed that: (i) the nitro- 
gen and carbon yields from low- and intermediate-mass 
stars are very uncertain and (ii) the carbon yields from 
stars with m > 40 Mq have been multiplied by a factor of 
three in order to mimic results from recent models taking 
mass loss and stellar rotation into account (see Chiappini 
et al. 2003a, b and references therein) and both of these 
mechanisms are still far from being fully understood. In 
conclusion, it is likely that the uncertainties in the stellar 
nucleosynthesis are the most important sources of errors as 
far as both carbon and nitrogen evolution are concerned. 
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The effect of extending the lognormal form of the IMF 
derived by Chabrier (2003) for m < 1 M© to the BD 
domain is also studied. BDs are low-mass, long-living ob- 
jects, that act simply as sinks of matter from the point 
of view of galactic chemical evolution. They were born at 
every time of Galaxy's evolution and never restored mat- 
ter into the ISM. According to what reported in Table 2, 
only a small fraction (< 10%) of the total mass of a stel- 
lar generation should be locked up in BDs. Therefore, 
only small effects are expected to appear when extend- 
ing the IMF to the BD domain. This is what is actu- 
ally found (compare model C03a to model COSb results 
in Table 5). Nevertheless, considering the existence of a 
substellar mass regime reduces the amount of matter pro- 
cessed through nuclear burning in stars during the whole 
Galactic lifetime. This results in lower abundances in the 
gaseous matter at the time of Sun's formation, thus ob- 
taining a better agreement with the observations. The 
current stellar density we obtain in the solar neighbour- 
hood in the two cases is fairly similar, being Sstars(iGai) — 
34 Mq pc~^ if BDs are not taken into account and 
Sstars(iGai) — 29 Mq pc~^ if they are. Both values are 
within the errors associated with the observational esti- 
mate, though the higher value is more likely. 

Figs. and El show the errors associated with chemi- 
cal evolution model predictions for several quantities, ob- 
tained by changing the prescriptions for the stellar IMF 
and keeping the same stellar lifetimes (Maeder & Meynet 
1989). Model predictions on the evolution of ^He/H and 
^He/''He as functions of time (Fig. upper and lower 
panel, respectively) and [0/Fe] and [S/Fe] as functions 
of [Fe/H] (Fig. El upper and lower panel, respectively) in 
the solar neighbourhood are shown, as well as the cor- 
responding data. Model predictions are not normalized 
to the predicted solar values, as instead is often done in 
chemical evolution studies. This is to better appreciate the 
differences obtained with the various IMFs in the last 4.5 
Gyr of evolution, an effect which is not evident when the 
solar normalization is applied. The elements displayed in 
Figs. and El are not chosen randomly. They are repre- 
sentative of stellar progenitors belonging to different ini- 
tial mass ranges: (i) ^He is produced by stars belonging 
to the low-mass range, 1-2 Mq; (ii) ^He comes from the 
whole stellar mass range; (iii) oxygen is almost entirely 
produced on short time scales by stars with m > 10 Mq; 
(iv) sulphur is representative of elements produced partly 
by SNell and partly by SNela; finally, (v) Fe is mostly 
produced by Type la SNe on long timescales. 

Scalo's (1998) IMF, having the highest fractions of 
stars with m > 10 M© and 1 < tu/Mq < 2 (Fig. HI), 
predicts a too large O /Fe ratio during the whole Galactic 
lifetime and an increase of ^He/H and ^He/'^He from the 
time of Sun's formation up to now that hardly fits the 
data. This is because the stars in the range 1-2 M© are 
net ■^He producers, even when ~90% of them are assumed 
to experience the cool bottom processes which strongly 
deplete their "^He yields (see next section for more details 
and references), whereas high- mass stars produce almost 
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Fig. 7. Type la SN rates obtained with different assump- 
tions on the stellar lifetimes (Maeder & Meynet 1989 ~ 
thin solid line; Tinsley 1980 - thick dotted line; Schaller 
et al. 1992 - thick solid line; Kodama 1997 ~ thick dashed 
line). We also show the effect of changing the fraction of 
mass entering the formation of Type la SN progenitors. 
The four curves lying in the lower part of the diagram 
have all been computed with A = 0.05, while the value of 
A for the curves lying in the upper part of the diagram 
has been changed so as to produce the same present-day 
_Ria in the disk, independently of the adopted stellar life- 
time prescriptions. The Type la SN rate observed in the 
Galaxy at the present time is also shown [vertical bar, for 
h = 0.7; Cappellaro et al. 1997). 



all the oxygen. On the contrary, Salpeter, Tinsley and 
Scalo (1986) IMFs predict enrichment histories for ^He/H 
and ^He/^He which better agree with the data, because of 
the lower mass fraction falling in the 1-2 M© stellar mass 
range. The Kroupa et al. and Chabrier IMFs, with x = 1.7 
for m > 1 Mq, produce an intermediate behaviour. They 
also guarantee the best fit to the [0/Fe] vs [Fe/H] data. 
The (small) differences in the S/Fe ratio in the disc are 
mostly related to the (slightly) different behaviour of the 
SNla rate in the past predicted when adopting different 
IMFs. 

5.2. Changing the stellar lifetimes 

In this section we keep fixed the prescriptions on the stellar 
IMF, while varying those on the stellar lifetimes, according 
to what has been reviewed in Section 3. We choose to 
adopt Scalo (1986) IMF for the sake of comparison with 
previous results published in a series of papers dealing 
with different aspects of the Milky Way evolution (e.g. 



10 



D. Romano et al.: Quantifying the uncertainties of diemical evoiution studies 



Table 6. Abundances in the PSC as predicted by the model at t = 9.2 Gyr for different assumptions on the stellar 
lifetimes. The abundances are by number in log£(X), except for helium (F) and global metallicity (Z), for which the 
abundances by mass are reported. Model predictions are compared to observed photospheric solar abundances, unless 
otherwise stated. A value oi A = 0.05 is assumed for all models (see text). 





Y 


C 


N 


O 


No 


Mg 


Si 


S 


Fc 


Z 


T80 


0.261 


8.39 


8.13 


8.90 


7.97 


7.23 


7.57 


7.18 


7.57 


0.017 


M89 


0.260 


8.36 


8.04 


8.77 


7.84 


7.14 


7.52 


7.11 


7.51 


0.014 


S92 


0.260 


8.36 


8.06 


8.80 


7.89 


7.17 


7.55 


7.10 


7.54 


0.015 


K97 


0.260 


8.35 


8.08 


8.81 


7.89 


7.17 


7.51 


7.11 


7.56 


0.015 


O bscrvations 


GS98 


0.275 ± 0.01 


8.52 ± 0.06 


7.92 ± 0.06 


8.83 ± 0.06 


8.08 ± 0.06" 


7.58 ± 0.05 


7.55 ± 0.05 


7.33 ± 0.11 


7.50 ± 0.05 


0.017 


HOI 




8.59 ± 0.11 


7.93 ± 0.11 


8.74 ± 0.08 




7.54 ± 0.06 


7.54 ± 0.05 




7.45 ± 0.08 




APOl 








8.69 ± 0.05 














AP02 




8.39 ± 0.04 


















A04 








8.66 ± 0.05 


7.84 ± 0.06 










0.0126 



Coronal data. 



T80 - Tinslcy (1980): M89 - Macdcr & Mcynct (1989); S92 - Schallor ct al. (1992): K97 - Kodama (1997); GS98 - Grcvcsso & Sauval (1998); HOI 
- Holwcgor (2001); APOl - AUcndc Pricto ct al. (2001); AP02 - Allcndc Pricto ct al. (2002); A04 - Asplund ct al. (2004). 



Chiappini et al. 1997; Romano et al. 2003, and references 
therein). 

In the following, we compare results of chemical evolu- 
tion models for the solar vicinity adopting stellar lifetimes 
from: (i) Tinsley (1980); (ii) Maeder & Meynet (1989)^; 
(iii) Schaller et al. (1992)^; (iv) Kodama (1997)'*. 

In Chiappini et al. (1997) and subsequent work by 
that group, the prescriptions of Maeder & Meynet (1989) 
were adopted. According to those authors, stars with 1.3 
< iti/Mq < 3 are characterized by lifetimes longer than 
in any other study (see Fig. |2J). This causes the Type la 
SN rate to increase smoothly during almost the whole 
disk evolution, with a gentle decline starting only a cou- 
ple of Gyrs ago (thin solid line in Fig. [T)). On the con- 
trary, models with different prescriptions on the stellar 
lifetimes all produce a well-defined peak in the rate, fol- 
lowed by a steep decline afterward (thick lines in Fig.EJ. 
If one adopts Tinsley's (1980), Schaller et al.'s (1992) or 
Kodama's (1997) stellar lifetimes, keeping A - i.e. the frac- 
tion of mass belonging to binary systems which will give 
rise to Type la SN explosions ~ constant and equal to 
0.05 results in a little bit more iron produced until the 
time of Sun's formation (see also Table 6, column 10), 
while SNela are supposed to be less numerous in the disk 
at the present time. However, their number still agrees 
with that inferred from the observations (vertical bar in 
Fig. Cj). To rise the present-day SNIa rate to the value 
expected when adopting Maeder & Meynet (1989) stellar 
lifetimes, it is necessary to increase the fraction of mass 
belonging to SNIa progenitors. In particular, when adopt- 
ing Schaller et al. (1992) or Kodama (1997) r^s, A = 0.08 
is required, while a slightly higher value, A = 0.085, must 
be assumed if Tinsley's (1980) stellar lifetimes are adopted 

^ In the m < 1.3 A/0, m > 60 Mq mass ranges we use the ex- 
trapolations adopted by Matteucci and coworkers (Matteucci 
& Frangois 1989; Chiappini et al. 1997). 

^ We use the polynomial fits of Gibson (1997). 

* His prescriptions are equivalent to those given in Padovani 
& Matteucci (1993). 



{upper curves in Fig. CJ. Increasing the efhciency of star 
formation while leaving the A parameter unchanged would 
also lead to an increase of the SNIa rate, but it would 
also cause major problems in reproducing other impor- 
tant constraints available for the immediate solar vicinity. 
From Fig. [7| and Table 6, column 10, we conclude that, 
for a Scalo (1986) IMF, a small A value of the order of 
0.05 or even less should be preferred, because it produces 
theoretical predictions matching remarkably well the solar 
iron abundance and the present-day SNIa rate data at the 
same time. 

The G-dwarf metallicity distributions obtained with 
all the above-mentioned prescriptions for Tm are displayed 
in Figs. ^ (upper left panel) and |H1 Again, the theoret- 
ical distributions are convolved with a Gaussian to ac- 
count for both the observational and the intrinsic scat- 
ter. The adopted total dispersion in [Fe/H] is tr = 0.15 
dex (Arimoto, Yoshii & Takahara 1992; Kotoneva et al. 
2002). Since most of the iron in the solar neighbourhood 
comes from Type la SNe, the shape of the distribution is 
expected to change according to changes affecting Type 
la SN progenitors. Indeed, adopting Tinsley (1980) or 
Kodama (1997) stellar lifetimes (Fig.|Sl left and right pan- 
els, respectively) results in a more pronounced peak, in- 
dependently of the value of A (which is given in the upper 
left corner of each panel, third row). 

It is worth noticing that, independently of the choice 
of the stellar lifetimes, we always get a strikingly good 
fit to the position of the peak and to the high-metallicity 
tail of the distribution. We conclude that: (i) the main 
parameter driving the location of the peak and the shape 
of the distribution is the adopted time scale for thin-disk 
formation, which was already well known (Chiappini et al. 
1997); (ii) the adopted stellar lifetimes affect the theoreti- 
cal G-dwarf distribution as well, though through a second 
order effect: they mostly act on the width and height of the 
distribution. However, convolution with a Gaussian which 
accounts for errors makes different distributions, obtained 
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with different prescriptions on tlic stellar lifetimes, look 
pretty much the same. 

In Table 6 we report the abundances in the Protosolar 
Cloud (PSC) predicted under different assumptions on the 
stellar lifetimes. Notice that for all the models in Table 6 
a value of 0.05 for A is adopted. A larger A value would 
overestimate the iron content of the Sun. It can be seen 
that the global metallicity, Z, is well reproduced by all the 
models and so are the CNO abundances, with a possible 
exception for oxygen, whose predicted solar abundance is a 
bit higher than observed. On the other hand, the helium 
abundance turns out to be too low. Increasing the star 
formation efficiency would provide a higher Y abundance, 
but the oxygen and global metallicity would also increase, 
breaking the agreement with the observations. This clearly 
shows that some revision on the helium stellar yields is 
necessary (see Meynet & Maeder 2002; Chiappini et al. 
2003b). 

Let us now comment on specific trends predicted for a 
handful of important species: (i) ^He and (ii) ^Li, coming 
mostly from low-mass stellar progenitors; (iii) oxygen, a 
typical massive star product; and (iv) sulphur, with both 
a SNII and a SNIa component. Major differences are ex- 
pected in the model predictions for '^He and ^Li, because of 
the large differences characterizing different r„j prescrip- 
tions for their stellar progenitors. 

To deal with '^He, one must keep in mind that it is 
mostly produced by low-mass - hence long-living - stars. 
Standard stellar evolutionary theory predicts a large pro- 
duction of ^He from low-mass stars leading to severe in- 
consistency between the observed '^He abundances and 
those predicted by chemical evolution models assuming 
these standard yields (Rood, Steigman & Tinsley 1976; 
GaUi et al. 1995; Olive et al. 1995; Dearborn et al. 1996; 
Prantzos 1996). A problem that has been now superseeded 
by the inclusion of rotational mixing in stellar models (e.g. 
Charbonnel 1995; Sackmann & Boothroyd 1999). Indeed, 
extra mixing makes it possible to reconcile predictions 
from Galactic chemical evolution models with observa- 
tions, as long as '^He is destroyed in a large enough fraction 
(~90%) of low-mass stars (e.g. Galli et al. 1997; Chiappini 
et al. 2002, and references therein). Fig. El shows the ^He 
evolution in the solar neighbourhood as predicted under 
different assumptions on the stellar lifetimes: the thin solid 
lines are for Maeder & Meynet (1989); the thick solid 
lines are for Schaller et al. (1992); the thick dotted lines 
for Tinsley (1980) and the thick dashed ones for Kodama 
(1997). The prescriptions on ^He synthesis are those from 
Dearborn et al. (1996) and Sackmann & Boothroyd (1999) 
for stars without and with extra-mixing. These "^He yields 
were recently adopted also by Chiappini et al. (2002) and 
Romano et al. (2003) and take extra-mixing in 93% of low- 
mass stars into account. The lower '^He production for t < 
11-12 Gyr predicted assuming Maeder & Meynet (1989) 
stellar lifetimes is due to the longer lifetimes in the mass 
range 1-2 Mq . The subsequent steep rise is mostly due to 
the late contribution from stars in the 0.6-0.9 M© mass 
range, which die if adopting the Matteucci & Frangois 
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Fig. 10. Temporal evolution of lithium in the solar neigh- 
bourhood. Different curves refer to model predictions 
obtained with different recipes for the stellar lifetimes. 
Thin solid line: Maeder & Meynet (1989); thick solid 
line: Schaller et al. (1992); thick dotted line: Tinsley 
(1980); thick dashed line: Kodama (1997). The meteoritic 
(Nichiporuk & Moore 1974; pentagon) and the depletion 
corrected value for the local ISM (Knauth et al. 2003; at 
1- and 2-tT - thick and thin vertical bars, respectively) are 
shown as well. Notice that the a value (where a is the 
fraction of white dwarfs entering the formation of nova 
systems) is adjusted so as to reproduce the nova outburst 
rate currently observed in the Galaxy. In particular, in 
the case of Tinsley, Schaller et al. and Kodama stellar 
lifetimes a ~ 0.02, while a ^ 0.01 if Maeder & Meynet 
stellar lifetimes are adopted instead (see text for details). 

(1989) extrapolation of Maeder & Meynet's stellar life- 
times in the very low stellar mass range. Conversely, these 
stars never die if the Schaller et al., Tinsley or Kodama 
stellar lifetimes are adopted instead (Fig. |21 see also the 
discussion in Romano et al. 2003). If one trusts the local 
value of the helium isotopic ratio as measured with the 
COLLISA experiment on board the space station MIR 
(Salerno et al. 2003 - Fig. 13 lower panel, vertical bar on 
the right at t = 13.7 Gyr), longer lifetimes for stars with 
m < 1 Mq should be preferred. In fact, this measure- 
ment supports the hypothesis that negligible changes of 
the abundance of '^He occurred in the Galaxy during the 
last 4.5 Gyr. 

Lithium is another element originating mostly from 
low-mass stars (Romano et al. 2001). According to our 
models, its meteoritic abundance is due for ^40% to low- 
mass (m = 1-2 Mq) stars on the red giant branch (RGB) 
and for ^10% to novae, binary systems consisting of a 
white dwarf plus a main sequence companion (Romano et 
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al. 2001, 2003; but see also Casuso & Beckman 2000 and 
Travaglio et al. 2001 for different views). The predicted 
steep rise of the lithium abundance in the ISM at late 
times is dictated by the long time scales for lithium pro- 
duction from RGB stars and novae when adopting Maeder 
& Meynet stellar lifetimes. Here we analyse what changes 
are introduced by adopting different stellar lifetimes. In 
Fig. 1101 we compare model predictions for the temporal 
evolution of ''Li in the solar vicinity obtained by adopting 
Maeder & Meynet (1989) stellar lifetimes (thin solid line) 
to what we get if adopting Schaller et al. (1992 - thick 
solid line), Tinsley (1980 - thick dotted line) or Kodama 
(1997- thick dashed line) prescriptions. The Schaller et al. 
stellar lifetimes result in a more gentle rise during the last 
3 Gyrs, while almost no evolution is expected during the 
same time interval with Tinsley (1980) or Kodama (1997) 
prescriptions. These considerations read interesting in the 
light of recent claims of a null ^Li evolution in the ISM 
during the last several Gyrs inferred from cluster and field 
star lithium data (Lambert & Reddy 2004). However, in 
our opinion these data do not rule out a scenario of late 
lithium pollution from low-mass stars, if the uncertainties 
in both the models and the data are properly taken into 
account. 

In our model, nearly 10% of lithium in meteorites is 
produced during nova outbursts. These are periodical ex- 
plosions, that do not destroy the parent system. When 
dealing with nova systems, one must introduce in the 
model a free parameter, a, describing the fraction of white 
dwarfs which enters the formation of nova systems (simi- 
larly to what is done for SNela; see D'Antona & Matteucci 
1991; Romano et al. 1999; Matteucci et al. 2003). Fig.Hni 
shows model results obtained by changing the a value so 
as to obtain the same theoretical nova outburst rate in 
the Galaxy at the present time, whatever the T„i choice. 
A value of Rnovaitoai) — 20 yr~^ is found with a ^ 0.01 in 
the case of Maeder & Meynet stellar lifetimes and a ^ 0.02 
in the remaining cases, to be compared with = 20- 

30 yr^i (Shaffer 1997). When using Tinsley's, Schaller 
et al.'s or Kodama's stellar lifetimes, a ^Li production 
from novae higher than expected when using Maeder & 
Meynet 's stellar lifetimes is obtained during the whole 
Galactic lifetime, owing to the higher value assumed for 
a. Lowering the ^Li yields from red giants and/or novae 
can make the theoretical predictions in better agreement 
with the observations. This goes in the right direction. In 
fact, in order to reproduce the observed meteoritic lithium 
data, in previous work we had to require that all low-mass 
red giants are displaying the largest observed atmospheric 
lithium enrichment near the tip of the RGB and couple 
this with the most efficient mass loss still compatible with 
observations (Romano et al. 2001). Adopting stellar life- 
times prescritions different from Maeder & Meynet (1989) 
helps us to alleviate such an extreme, ad hoc scenario. 

Finally, as far as the [0/Fe] and [S/Fe] ratios are con- 
cerned, it is worth emphasizing that only small differences 
are found among models adopting different stellar life- 
times. In the case of Tinsley's (1980) prescriptions a flatter 



behaviour is predicted for [0/Fe] vs [Fe/H] at high metal- 
licities, at variance with observations (Bensby, Feltzing, 
& Lundstrom 2004). However, no firm conclusions can be 
drawn on this point alone. 

6. Final remarks and conclusions 

Together with the time modulation of the SFR, the IMF 
dictates the evolution and fate of galaxies. Nonluminous 
BDs and the lowest-mass stars lock up an increasing frac- 
tion of the baryonic mass of galaxies over the cosmologi- 
cal time. Short-lived massive stars and intermediate-mass 
stars belonging to binary systems ending up as Type la 
SNe heat the ISM, eventually determining the occurrence 
of galactic outflows. It is therefore of much importance 
to quantify the effect of changing the relative numbers of 
stars in different mass ranges at different times on galactic 
chemical evolution model predictions. 

In this work we deal with our own Galaxy. We as- 
certain the range of variations that affect several model 
predictions when accounting for uncertainties on both the 
stellar IMF and the stellar lifetimes. First, we show results 
of chemical evolution models for the Milky Way computed 
with different assumptions on the stellar IMF. Then, we 
investigate the effects of changing the prescriptions on 
the stellar lifetimes. 'Theoretical errors' are associated to 
model predictions due to uncertainties in the IMF and/or 
stellar lifetimes. We summarize our main conclusions as 
follows: 

Among all the studied IMFs, the Salpeter (1955) and 
Scalo (1998) ones are those in less agreement with the 
data. The Scalo (1986), Kroupaet al. (1993) and Chabrier 
(2003) IMFs all guarantee the best flts to several impor- 
tant observed properties of the solar vicinity. 

Different stellar lifetime prescriptions differ mostly in 
the low and very low stellar mass domain. Therefore, it is 
not surprising that models adopting different prescriptions 
for the stellar lifetimes differ mostly in the predicted evo- 
lution for those species originating mostly from low-mass 
stars. We analyse the evolution of '^He and ^Li and show 
that it is better reproduced if adopting long lifetimes for 
the low stellar mass range. 

Oxygen abundance data recently derived for F and G 
dwarf stars in the solar neighbourhood indicate that the 
[0/Fe] trend at super-solar [Fe/H] continues downward, 
which is in concordance with models of Galactic chemical 
evolution unless stellar lifetimes from Tinsley (1980) are 
adopted. Notice that her prescriptions are characterized 
by very short lifetimes for massive stars. 

We conclude that, given the uncertainties still associ- 
ated to current observations, the main conclusions reached 
by chemical evolution models for the Galaxy are left un- 
changed. However, it is clear that only further observa- 
tions and studying galaxies of different morphological type 
will allow us to draw a flrmer picture and a better under- 
standing of the problem. 
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Fig. 4. Theoretical G-dwarf metallicity distribution as predicted by the model with Maeder & Meynet (1989) stellar 
lifetimes and different assumptions on the IMF (thick lines) compared to the data (histograms; solid: Rocha-Pinto & 
Maciel 1996; dotted: J0rgensen 2000). Theoretical [Fe/H] are normalized to the iron to hydrogen ratios predicted by 
the model at Sun's birth, which are displayed in the upper left corner of each panel. The first value refers to a common 
choice for the parameter regulating the fraction of stars giving rise to SNIa events, A = 0.05; the second value refers to 
the case in which the A parameter is adjusted so as to predict the same SNIa rate at the present time when changing 
the IMF (see text for details). The theoretical distributions corresponding to the latter case are shown as dashed 
areas. Remember that SNela are the major iron producers in the Milky Way. Notice that theoretical distributions are 
convolved with a Gaussian in order to account for both the observational and the intrinsic scatter. 
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Fig. 5. Upper panel: evolution of ^He/H in the solar neighbourhood. Different Unes refer to different IMF prescriptions: 
short-dashed line: Salpeter (1955); dotted line: Tinsley (1980); continuous line: Scalo (1986); long-dashed line: Kroupa 
et al. (1993); dot- short- dashed line: Scalo (1998); dot-long- dashed line: Chabrier (2003) with x = 1.7 for m > 1 Mq. 
Data (vertical bars at t — 9.2 and 13.7 Gyr) are from Geiss & Gloeckler (1998). Lower panel: evolution of ^He/*He in 
the solar neighbourhood. Models are labeled as in the upper panel. Data are from Geiss & Gloeckler (1998; at 1 and 
2a - thick and thin vertical bars, respectively, at t = 9.2 and 13.7 Gyr) and Salerno et al. (2003; 1 a-bar on the right 
at t = 13.7 Gyr). 
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Fig. 6. [0/Fe] (upper panel) and [S/Fe] (lower panel) ratios as functions of [Fe/H]. Models as in Fig.|Hl All the ratios 
are normalized to solar elemental abundances by Grevesse & Sauval (1998), except for oxygen, for which the solar 
value from Holweger (2001) is adopted. Data for oxygen are from Melendez & Barbuy (2002; stars: [0/Fe] values 
from infrared OH lines; full circles: mean [0/Fe] from [Of] line in bins of 0.2 dex in [Fe/H]; the size of the circles 
represents the number of stars in each metallicity bin). Data for sulphur are from Ryde & Lambert (2004; pentagons: 
measurements from observations of Sf hues at 9212.9, 9228.1 and 9237.5 A) and Chen et al. (2002; crosses). For these 
latter measurements, the typical error is also reported on the lower right corner of the panel. 
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Fig. 8. Theoretical G-dwarf metallicity distributions (thick lines) predicted by the model with different assumptions 
on the stellar lifetimes {left panels: Tinsley 1980; middle panels: Shaller et al. 1992; right panels: Kodama 1997) and 
on the fraction of mass entering the formation of Type la SN progenitors (upper panels: A = 0.05; lower panels: A = 
0.08 in case of Schaller et al.'s and Kodama's stellar lifetimes; A = 0.085 in case of Tinsley's stellar lifetimes). Scale's 
(1986) IMF is assumed in all cases. Theoretical distributions are compared to observational ones from Rocha-Pinto & 
Maciel (1996; solid histograms) and J0rgensen (2000; dotted histograms). For each model, the theoretical [Fe/H] ratios 
are normalized to the Fe/H ratio at Sun's birth predicted by the model itself, given in the upper left corner of each 
panel (middle row). This produces a shift of the distribution along the x axis, thus allowing for a better comparison of 
the theoretical shape with the observed one even if the model does not recover the actual iron abundance of the Sun. 
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Fig. 9. Same as Fig. but here different Hne types refer to different prescriptions on the stellar lifetimes. Thin solid 
lines: Maeder & Meynet (1989); thick solid lines: Schaller et al. (1992); thick dotted lines: Tinsley (1980); thick dashed 
lines: Kodama (1997). The Scalo (1986) IMF is assumed for all the models. 



